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Conjugated polymers are soft, one-dimensional conductors that admit complex interactions be- 
tween their polymeric, conformational degrees of freedom and their electronic ones. The presence of 
extended electronic states along their backbone allows for inter-chain electronic tunneling at points 
where these polymers make near passes. Using a combination of analytic modeling and Hartree-Fock 
numerical calculations, we study the localized electronic states that form due to such close encoun- 
ters between semiconducting conjugated polymers and explore how these states lead to chain-chain 
binding. We also study the interaction of these inter-chain binding sites with solitonic excitations 
on the chains. From these results and a modified Poland-Scheraga model, we determine the equi- 
librium structures of paired-chains formed by intermolecular electronic tunneling. We calculate the 
energetic ground state of such pairs and show the effective thermal persistence length of the paired 
chains can vary over an order of magnitude due to the intermolecular binding mechanism. 

PACS numbers: 82.35.Cd,36.20.Kd,72.80.Le 



I. INTRODUCTION 

Conjugated polymers [iHSl can be thought of as soft, 
low-dim.ensional conductors or semiconductors. Due to 
their promise in technological applications such as LEDs, 
solar cells Q-Ql , and even biosensors 0, Q , a great variety 
of such materials have been synthesized and studied. In- 
dependent of their specific chemical details, the existence 
of extended electronic states along their polymeric back- 
bone provides for their surprising (for a polymer) elec- 
tronic properties. We refer to these (semi-)conductors as 
"soft" since their polymeric backbones typically have a 
short (nanometer scale) thermal persistence length; the 
conduction path along the polymers in equilibrium solu- 
tion should be considered to be a tortuous random-walk. 

These molecules admit strong correlations between 
their fluctuating local configurational and electronic 
properties Q. In the extreme case, sharp (i.e. local- 
ized) bends in the polymeric backbone result in regions 
of poor electron transport at these so-called conjugation 
breaks. Conversely, the electronic degrees of freedom of 
the molecule influence its conformational degrees of free- 
dom. High electronic mobility along the thermally fluctu- 
ating backbone's random path allows for complexities in 
the interactions of a conjugated polymer with itself and 
with other such molecules. For example, in conducting 
chains the pressure of the electron gas trapped between 
nearby conjugation breaks (sharp bends) in the chain 
should force such sharp bends apart and thus enhance 
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the statistical weight of locally straighter backbone con- 
tours [l^. Thus, the electronic degrees of freedom con- 
tribute to the effective persistence length of the polymer. 
Moreover, looped conformations of the polymer backbone 
generating close contacts of parts of the polymer sepa- 
rated by a large arc length along its backbone may further 
perturb the electronic structure of the molecule by intro- 
ducing tunneling between distant sites along the chain. 
This sort of "bridge conduction" [ll| leads to localized 
attractive interactions of the polymer at such crossing 
points where the chain makes a near pass to either it- 
self or another. Based on this electronically mediated 
intra-molecular attraction, one may imagine that conju- 
gated polymers can self-aggregate into compact globules 
in solution [ll|. At higher polymeric concentrations, one 
would expect to observe interchain binding, the subject 
of this article. 

More recently, detailed numerical calculations have 
been preformed that determined the magnitude of these 
inter-chain tunneling matrix elements for specific con- 
jugated polymers [i2|. Using these results and a simple 
tight-binding model for metallic [1^ and/or semiconduct- 
ing [Tjl chains, we showed that, for the case of pairs of 
chains, intermolecular tunneling leads to attractive in- 
teractions for which the binding energy is ^ fcsT. The 
tunneling creates a pair of electronic bound states local- 
ized at the inter-chain crossing, leading to a total de- 
crease of the electronic contribution to the system's en- 
ergy on this scale. One may look at this attractive inter- 
action as an analogue to the traditional covalent chemi- 
cal bond, but at an energy scale two orders of magnitude 
weaker. We understand the weakness of the bond to re- 
sult from a combination of molecular geometry creating 
a larger than typical internuclear distance in the bond 



and the energetic cost of localizing the electrons from 
the extended states of the molecule in order to fill the 
new localized bound state. The magnitude of this in- 
termolecular, short-ranged interaction suggests that the 
statistical mechanics of intermolecular binding driven by 
this mechanism is rather subtle, since both chain config- 
urational entropy and chain bending energies contribute 
to the free energy of the system at the same energy scale. 

In the current article we study in more detail the effect 
of inter-chain tunneling sites on the electronic states of 
semiconducting chains by exploring analytically the in- 
teraction of pairs of such tunneling sites, the interaction 
of solitonic chain excitations with such tunneling sites, 
and determining the electronic ground state of pairs of 
chains bound uniformly along their length. We also use 
numerical calculations incorporating more details of the 
chemical structure of the molecules in question to study 
the dependence of the electronic binding interaction on 
the local geometry of the crossing point between the two 
chains. Finally, we use these results as input to a gener- 
alized Poland- Schcraga model describing the binding of 
two chains. Using this model, we determine the critical 
concentration for chain pairing and study how the con- 
formational statistics, as parameterized by the effective 
thermal persistence length of the paired chains evolves as 
a function of the strength of the binding interaction. 

To introduce the bandgap at the Fermi level and 
thus create the semiconducting state, we use the SSH 
model [l| for polyacetylene. This simple tight-binding 
system undergoes a Peierls instability [1^. While this 
model is based on the specific and uniquely simple model 
of polyacetylene, we believe that the results we obtain 
apply more broadly to semiconducting conjugated poly- 
mers. Indeed, similar tight-binding models have proven 
useful in predicting not only the properties of polyacety- 
lene, but have been found to be a useful model for 
more chemically complex conjugated polymers such as 
PPV [la: 113 ■ Thus, we present this model as a general 
one for the interactions between the electronic degrees of 
freedom and the conformational state of the polymer. 

Our principal results may be summarized as follows. 
We find (section III Ap that the magnitude of the single- 
site binding energy and the nature of the bound elec- 
tron wavefunction are both very similar for the bind- 
ing metallic or semiconducting chains. This similarity 
is due to the fact that the binding energy is associated 
with the creation of localized states having energies far 
from the Fermi energy, so that the density of states at 
the Fermi energy in the unperturbed system is not rele- 
vant to the inter-chain interaction energy. These results 
were discussed previously in Ref. [ij. We also exam- 
ine tunneling site interactions with each other (section 
IIIBp and with solitons that are necessarily present on 
sufficiently long dimerized chains f section ITV)) . We also 
compare the lowest-energy bound state configurations of 
semiconducting and metallic chains. Although these two 
systems seem similar at the level of individual tunneling 
sites, in section Hill we find that the ground states differ 



rather dramatically [1J| . We also discuss the implications 
of our numerical quantum chemical calculations for the 
stability of these ground states, determined by the SSH 
model. In section|V]we employ these results in the devel- 
opment of a modified Poland-Scheraga model to describe 
inter-chain binding and the conformational statistics of 
the result bound chain pairs. We conclude in section IVTl 
and in the Appendices we present the details of our cal- 
culations and numerical quantum chemical calculations 
supporting our binding mechanism. 



II. ISOLATED BINDING SITES 
A. A Single Crossing Point 

The binding between conjugated polymers is a generic 
property of the mixing of partially filled 7r-orbitals. We 
treat the electronic component of the interaction using a 
simple tight-binding model. To introduce a band gap we 
use the elegant SSH formalism where a Peierls instability 
opens a band gap through a distortion of the one dimen- 
sional lattice. One could introduce similar tight-binding 
models with a more complex unit cell appropriate for 
e.g. PPV [l^, but such calculations mask the underlying 
physics of the system in its most straight forward version. 
This simplification of the problem has been exploited pre- 
viously by Guo and collaborators, who showed that the 
SSH model can be used to study other polymers, such 
as PPV, with a suitable reinterpretation of the model 
parameters [13 . Of course, these chemically distinct sys- 
tems would produce quantitatively different binding in- 
teractions, but our primary findings, e.g. that there is 
an attractive interaction due to the creation of localized 
states above and below each band, should be generally 
insensitive to their chemical details. The SSH model also 
allows us to explore in an analytically tractable manner 
the interaction between inter-chain binding sites and lo- 
cal lattice distortions such as solitons. These results are 
also generalizablc to other conjugated polymer systems 
lacking this ground state degeneracy as superpositions 
of soliton states can be used to model the more general 
polaron problem in these systems |20| . 

The SSH Hamiltonian is given by 

HssH - -^t„+i(|£+l,a)(£,a| + |£,a)(£+l,a|) 
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Here and throughout this article we neglect the dynamics 
of phonons not associated with the soft phonon mode at 
the edge of Brillouin zone. We will not consider electron- 
phonon scattering, but we will allow for more complex 
static lattice distortions when computing the interaction 
of tunneling sites and solitons. Here \i,(j) is the tight- 
binding state at the ith site on the chain with spin <t, 
K is the effective spring constant of the a bonds, and 



ug is the displacement of the f th carbon from its equihb- 
rium position at £a. Introducing the Peierls instabihty 
for the system with one electron per site we generate 
static lattice distortions Ui = {—lYu (on the scale of 
u ~ 0.04 A) and shift the hopping matrix elements so 
that tej+i = to — a{ue+i ~ ue), doubling the unit cell 
of the lattice. Now the electronic states of the system 
satisfy the dispersion relation (see Fig. [T^) 
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which has a total bandwidth of 4<o and a bandgap cen- 
tered around E = oi width 8uoa = 2A [y. 

To study the interaction between two chains, we add 
a second chain with the same Hamiltonian as Eq. [1] and 
an interaction term of the form 
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-t'(|0,l,a)(0,2,a| + |0,2,a)(0,l,a|), (3) 



where we have defined the origin to lie at the interaction 
site, and we have labeled the chains as 1 and 2. The inter- 
chain hopping parameter t' implicitly includes the effects 
of the medium and the relative orientation of the chains. 
It is also a function of the distance between the linked 
tight-binding sites with a decay length on the order of a 
Bohr radius. Quantum chemical calculations on similar 
conjugated systems have shown that t' is on the order of 
0.1 eV [12|. We assume the chains have adopted a confor- 
mation where only one pair of sites is close enough for the 
inter-chain hopping to be significant. This configuration 
is illustrated in Fig. [TJd and c. 

The interaction part of the Hamiltonian given in Eq. [3] 
can be diagonalized using basis of states that are symmet- 
ric or anti-symmetric combinations of the single-chain 
states, 

K,s,a) = 2-'/\\£,l,a) + \i,2,a)) 

\£,a,a) = 2-1/2(1^, 1, a) -K, 2, a)) (4) 

In this basis the Hamiltonian of the system breaks up 
into two copies of the single chain Hamiltonian (Eq. [1]) 
that are identical except for the presence of a diagonal 
element in each sub-space of the Hamiltonian at the in- 
teraction site with the value —t'{t') for the symmetric 
(anti-symmetric) sub-space. 

The tunnel matrix element is now analogous to that 
of an impurity atom in a one dimensional crystal. The 
impurity potential gives rise to bound states by scatter- 
ing growing states into decaying states. In appendix |X] 
we use a transfer matrix technique to show that the in- 
teraction leads to the creation of four bound states as 
illustrated by the (red) dots in Fig. [T^. Two of these 
states, termed "ultraband" states [21|, appear above the 
conduction band and below the valance band with ener- 
gies 
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FIG. 1: (color online)(top) The spectrum of electronic states 
from Eq. [2]is plotted against wavenumber fc, shown as (black) 
solid lines, along with the four bound states at k/2a — 0, tt, 
shown as (red) filled circles, resulting from the single site 
interaction Eq. [S] (middle) Cartoon of two polyacetylene 
molecules with one tunneling site (green) where there is a 
significant overlap of the pz molecular orbitals. The carbon 
atoms are shown as (blue) spheres connected by (red) lines 
representing the sigma bonds, (bottom) Schematic of the 
Hamiltonian that describes this configuration. The double 
(=) and single (-) lines correspond to the matrix elements ii 
and i2 respectively. The dotted line indicates the inter-chain 
interaction (i'). 



while the other two appear in the gap and at the edges 
of the Brillouin zone with energies 
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If we assume that t' , A ^ tQ wc find that, to lowest or- 
der, the ultraband states have been shifted away from 
the band edge by i'^/4to which is unchanged from the 



metallic case [l3|. In contrast, the gap states are shifted 
by i'^ A/Sip. For half- filled chains, only the two bound 
states associated with the valence band will be filled. Be- 
cause the shift of the gap state is smaller than the shift 
of the ultraband state by 0(A/io) for typical parameter 
values, the result is a net lowering of the energy. 

The presence of the binding site will also perturb the 
energy of the extended electron states. However, as in the 
metallic case, it can be shown that the total contribution 
of these states is 0{N~^), where N is the total number 
of tight-binding sites on the polymer, and therefore can 
be neglected for long chains jl3| . 

In appendix IB] we present Hartree-Fock (HF) calcula- 
tions that verify the qualitative features of this binding 
mechanism. These calculations show that the binding 
energy is of order ksT and that it is sufficiently insen- 
sitive to the precise orientation of the polymers to allow 
for a variety of aggregate morphologies. 



B. Two Binding Sites 

Since conjugated polymers are semiflexible, it is pos- 
sible for them to bind in multiple places separated by 
arc lengths of unbound chains. Two such distant binding 
sites will experience an effective interaction due to a com- 
bination of the change in chain configurational entropy 
due to the binding constraint, the bending elasticity of 
the chains, and possibility the electrostatic repulsion of 
the backbones [24|. In addition, there may be a mod- 
ification of the binding energy of each binding site due 
to the interaction of the electronic wavefunctions associ- 
ated with each of them. We now study this latter effect 
by determining the shift in total bound state energy of 
two binding sites as a function of their separation along 
the chain. 

Previously, we demonstrated that the energy shift of 
an extended state in the presence of a single impurity 
could be expressed as a series in t' /IqN [13|. Because 
the symmetric and anti-symmetric states differ only in 
the sign of the impurity, the odd terms in this series 
will cancel. Therefore, after summing over the 0{N) 
extended states, the total energy contribution from these 
states must vanish as N~^. 

A second binding site introduces a new length d, the 
distance between the binding sites. We expect the per- 
turbation series for the extended states to contain terms 
proportional to d~^, suggesting that the contribution of 
the extended states may be non- vanishing for long chains. 
This is indeed the case. Taking these effects into account, 
we determine the interaction energy of two binding sites 
by direct numerical diagonalization of Eqs. lA3l the result 
is shown in Fig. [51 There is an attractive potential well 
for the two tunneling sites with a minimum at the sepa- 
ration of 2a. At separations greater than about 10 tight 
binding sites the binding energy is essentially constant, 
however, at all smaller separations there is a net reduc- 
tion in energy for even separations and a corresponding 



increase in energy for odd separations. 
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FIG. 2: The interaction potential of two tunneling sites as a 
function of the distance d between them measured in units 
of the unperturbed lattice a, as determined by the numerical 
diagonalization of chains 100 tight binding sites each. The 
tunneling matrix elements are given by t' = 0.25 eV. For these 
parameters, the energy difference between bound d/a — 2 and 
free d/a > 10 tunneling sites is approximately O.lfcsT. This 
weak attraction increases with t' . Entropic effects, however, 
play a comparable role in binding site interactions. 

This "even-odd" effect can be understood by looking at 
the states that receive the maximum perturbation from 
the binding potential. From first order perturbation the- 
ory, we expect that the energy shift of a given state to 
be proportional to the amplitude of that state at the two 
impurities. Therefore, the states perturbed the most by 
the two impurities separated by distance d will be those 
with a dominant wavenumber k that satisfies the relation 
kd = mr. However, the interaction potential induces a 
phase shift, (j) in the wavefunction. So, the condition for 
the maximally perturbed state is actually kd + (j) = nn. 
Because (f) is positive for attractive impurity potentials 
and negative for repulsive potentials (23| , the symmetric 
states will see the full effect of the impurity potential at 
smaller values of k than the anti-symmetric states. 

For the present case of half-filled electronic states, the 
highest occupied level has a wavenumber kf = 7r/2a. For 
even d/a we find that kfd = mr so the states with large 
negative shifts are filled but the corresponding states with 
large positive shifts are unoccupied. However, if d/a is 
odd the states with large positive shifts are also filled 
resulting in an increase in energy relative to the two iso- 
lated binding sites. The interaction between two binding 
sites is discussed in more detail in Appendix [C] Given 
these results, we expect that binding sites should have 
a weak tendency to cluster so that chains develop finite 
bound sections in thermal equilibrium. Any conclusions 
regarding bound-state chain configurations will have to 
be postponed until we consider the chain configurational 
entropy of these bound states. We turn to this problem 
in scction|Vl Before we do so, we consider the problem of 
multiple inter-chain binding sites (instead of just two) as 



might occur when the chains adopt a parallel alignment. 



III. BINDING/NON-BINDING SYMMETRY 

If two chains are aligned in parallel, it is possible to 
have binding events at every site along the chain as illus- 
trated in Fig. [21 The interaction Hamiltonian is then 
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t'> (|£,l,a)(£,2,a| + |£,2,a)(^,l,<7|). (7) 






This interaction can also be diagonalized using the 
change of basis shown in Eq. |4l We find that both the 
valence and conduction bands are split with the chain- 
symmetrized states dropping in energy by t', and the 
chain-anti-symmetrized ones rising by t'; see Fig. [31 For 
conjugated polymers interacting in solution, however, t' 
is considerably smaller than A. Thus, the symmetric and 
anti-symmetric bands have equal occupancy, and there is 
no net change in energy and no binding. 

At first it seems counter-intuitive that the chains can 
bind at a single site, but have no attraction with many 
binding sites. To understand this result we now show 
that there is a symmetry between chains having spatially 
uniform densities of binding sites n and 1 — n. The lack of 
binding at maximal tunneling site density n = 1 \s thus 
understandable since this state is energetically equivalent 
to no binding sites at all: n = 0. This symmetry results 
in a maximally bound state with interchain tunneling at 
every other site, i.e. n = 1/2 [l4| . 

To observe this proposed symmetry, consider first the 
Hamiltonian of chains bound at every site: Hps = 
Hq + H'j where ifo is the symmetrized/anti-symmetrized 
non-interacting Hamiltonian, and H'j is the symmetrized 
version of Eq. [71 As shown in Fig. [J]:, this Hamiltonian 
is identical to Hq apart from equal and opposite constant 
shifts in the symmetric and anti-symmetric bands. If we 
now remove the overlap at a single site, we have effec- 
tively introduced an (non-binding) impurity site having 
a potential of strength +i'(— i') in the symmetric (anti- 
symmetric) band. Because of the symmetry t' — > —t' in 
Eqs. [51 and [HI the removal of an interaction at a single 
site results in the same net energy change as adding a 
single interaction site to the non-interacting chains. This 
argument can be extended to any set of interacting sites. 
Allowing a set [a] of sites to interact on the unbound 
chains will change the energy of the symmetric band 
by the same amount as the energy change of the anti- 
symmetric band upon the removal of interactions from 
the same set {a\ from the fully bound chains. Similarly, 
the anti-symmetric band on the non-interacting chains 
changes by the same energy as the symmetric band on 
the fully bound chains upon the addition/subtraction of 
interactions at {a}. 

The electronic ground state of the pair of parallel 
chains, as a consequence of this symmetry, requires tun- 
neling at every-other site. The all-trans configuration 
of the carbon backbone of the polyacetylene molecule as 
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FIG. 3: (color online) (top) Cartoon of the paired chains hav- 
ing tunneling sites at every tight binding site. The color 
scheme is identical to that of Fig.[l] (middle) Schematic repre- 
sentation of the Hamiltonian describing same set of tunneling 
sites, (bottom) The electronic band structure for the paired 
chains with tunneling at every site. 



shown in Fig.^lallows for precisely the required geometry. 
It is interesting to note that this tunneling site pattern 
effectively widens the Peierls gap. 

More complex conjugated polymers may not be able 
to achieve this maximally bound state due to e.g. steric 
interactions of side chains. In spite of the symmetry be- 
tween pairs of chains having binding density n and 1 — n, 
we do not expect to observe cases where n > 1/2 as 
a reduction of the number of inter-chain tunneling sites 
will simultaneously decrease the electronic contribution 
to the energy of the system and increase the its eonfigu- 
rational entropy. 



given by the solutions of cubic polynomial 



= E^ - E{Ul + A^ + t'2) ± WtoA. 



(8) 





FIG. 4: (color online) (top) Cartoon of the paired chains hav- 
ing tunneling sites at every other tight binding site, using the 
color scheme of Fig.[l] This is the tunneling configuration that 
results in the lowest electronic energy of the pair by increas- 
ing the Peierls gap. (bottom) Schematic of the Hamiltonian 
describing this configuration. 



IV. SOLITONS 

The dimerization pattern of the tight-binding sites in 
the SSH model breaks a two- fold symmetry of the system. 
Reading from the left to the right the slightly longer bond 
(— ) can either precede or follow the slightly shorter one 
{—). In thermal equilibrium this one-dimensional system 
cannot maintain long-range order by keeping only one of 
these two dimerization patterns. One expects to find do- 
main walls where the dimerization pattern switches from 
= — = to o — = — o. At the domain wall, or soli- 
ton, the pattern of bond lengths is necessarily distorted 
over some finite distance. These solitonic structures have 
been studied in detail by Su and collaborators [l|: They 
found that there are new localized electronic states as- 
sociated with these doiTiain walls. It remains to be seen 
how these localized states associated with the solitons are 
affected by inter-chain tunneling. 

For computational simplicity, we first consider the 
symmetric case in which the two chains each have a sin- 
gle soliton located at the tunneling site. This allows 
us to separate the binding and anti-binding states us- 
ing Eq. H] Furthermore, we simplify the structure of the 
soliton m to that a single-site domain wall. Using our 
notational short hand, this soliton can be represented as 
= — = — • — = — = where the filled circle 
represents both the center of the soliton and the location 
of the binding site. In Appendix |X] we show that the 
bound states associated with this scattering center are 



In the limit that E ^ t' , Awe find roots corresponding to 
ultraband bound states at energies ±(4t§ + A^ -I- 1'^)^/^, 
which, except for the A^ term, are identical to the re- 
sult for metallic chains [13. We also find another root 
of Eq. [5] in the gap. This mid-gap state has an energy 
of 4t'toA/{4:tl -f A^ -f t'2). The presence of two solitons 
has minimal effect on the ultraband bound states, while 
the mid-gap states associated with the solitons are split 
by 0{t'A/to). For half-filled chains the lower gap state is 
filled for both spin states, while the corresponding upper 
state remains empty in the electronic ground state of the 
system. Although solitonic perturbation is linear in t', as 
opposed to the ultraband energy shifts which are ©(t'^), 
the net energy change is comparable to the binding en- 
ergy from ultraband states due to the additional factor of 
A/to- The net result is an 0(1) enhancement to the bind- 
ing energy due to the co-localization of the solitons at the 
inter-chain tunneling site. It is important to note that 
the binding energy enhancement associated with solitonic 
co-localization occurs only for uncharged solitons, since 
charged ones contribute either empty or doubly occupied 
mid-gap states. The splitting of the filled mid-gap states 
by the tunneling matrix element results in no net en- 
ergy change of the system. Now, the energy required to 
form a soliton is approximately 0.42eV [l| suggesting the 
thermal equilibrium density of uncharged solitons is typ- 
ically small, but they can also be trapped on long chains 
by the process of cis- to trans-isomerization; addition- 
ally a single uncharged soliton will spontaneously form 
on chains with an odd number of sites |24l . [25| . Thus, 
soliton co-localization can play a role in the strengthen- 
ing of isolated inter-chain binding sites, but presumably 
is less relevant as the number of binding sites increases. 

We studied the co-localization of uncharged solitons 
with the tunneling site numerically in order to explore 
a more physical extended soliton. Using a tanh[s/(As)] 
profile for the static displacement field of the soliton, we 
verified the 0(1) enhancement of the binding energy. The 
calculation was performed on two chains of 199 sites each 
with a soliton centered at site 101. We also varied the 
width of the soliton As to ascertain if the presence of 
the tunneling site altered the minimum energy structure 
of the soliton as determined by Su and collaborators [l|] . 
We found that such an effect was negligible. 

More generally, we determined the change in electronic 
ground state energy of two chains forming one tunnel 
junction and each having one soliton whose center is lo- 
cated Axi sites from the tunneling site. Here the index 
i labels the chains. These results are shown in Fig. \5\ 
The principal result is that both solitons are strongly at- 
tracted to the the tunneling site, but the strength of the 
attractive potential of a soliton on one chain is signifi- 
cantly enhanced when the soliton on the other chain is 
already localized at the tunneling site. This effect results 
from the mixing of the mid-gap states on each chain due 
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FIG. 5: a) (color online) The single tunneling site binding en- 
ergy of two chains (labeled (1) and (2)), each having a soliton 
centered at a distance Axi, Ax2 from the tunneling site. The 
minimum energy (maximum binding) occurs when both soli- 
tons are localized at the tunneling site: Axi = Aa;2 = 0. Due 
to the "even-odd" effect mentioned in the text, we represent 
the results as two surfaces, (top) The binding energy for cases 
where the distances Aa;i , Aa;2 represent even numbers of in- 
teratomic distances, (middle) The corrugated energy surface 
showing the binding energy for separations where at least one 
of Aa:! or Ax2 are odd. Note the change in energy scale from 
above, (bottom) A schematic representation of the centers 
of solitons at distances Aa;i and Ax2 from the tunneling site 
(dotted line). The figures were calculated using A*' = 199, 
4fo = 10 eV, uo = 0.04 A, K ^ 21 eV/A^ a = 4.1 eV/A, 
t' = 0.25 eV, and As = 7a. 



to the inter-chain tunneling matrix element and leads to 
a further reduction of the energy of the filled mid-gap 
states. We also note the large amplitude oscillation of 
the solitonic interaction energy with distance from the 
tunneling site. This is simply due to the fact that the 
wavcfunction of the mid-gap states has nodes on alter- 
nating tight-binding sites. 



V. POLYMER DEGREES OF FREEDOM 

For polymers in solution, the morphology of their ag- 
gregates is determined by the relative strength of the in- 
tcrmolccular binding mechanism and the conformational 
degrees of freedom of the polymers. These parameters 
will depend strongly on variables like electrostatic screen- 
ing and solvent mediation of the interchain hopping pa- 
rameter. Since these are expected to vary greatly system- 
by-system, we present a general calculation describing 
the onset of aggregation. 

In sufficiently dilute solution, we expect the formation 
of stable two-polymer aggregates. These may take the 
form of loosely-bound, braid-like structures in which the 
polymers form non-interacting loops between rare bind- 
ing sites, or zipper-like with the polymers interacting at 
many consecutive sites j22| . Since the single-site binding 
mechanism discussed in section III Al is qualitatively un- 
changed from the doped case, the braid structures will 
also be unchanged from conducting polymers discussed 
earlier [l^ . The formation of the tightly-bound, so-called 
"zipper structures," however, differs between the con- 
ducting and semiconducting polymers. As shown in that 
previous work on the doped, metallic system, the length 
of the tightly bound regions is limited by the doping level 
of the polymers. In the case of undoped, semiconduct- 
ing polymers, the binding energy due to long, tightly- 
bound regions grows linearly with their length, provided 
the binding sites occur at every other site as shown in 
Fig. 21 For a pair of bound polymers interacting via a 
series of well-separated tightly-bound regions, the total 
density of states will be equal to the sum of the den- 
sity of states for independent segments of tightly bound 
and free chains weighted by the respective lengths. For 
this situation (with binding sites at every other site) the 
dispersion relations for the conduction bands are given 
by 
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the valance bands have equal magnitude and opposite 
signs. We find that the binding energy per site in the 
tightly bound segments is then given by 
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where = ka and additional factors of two have been in- 
cluded to account for spin, positive and negative values of 
k, and the fact that two chains are present. This binding 
energy is the product of the same dimcnsionful prefactor 
t'^/2to, found for the single site binding energy, and a 



dimcnsionlcss integral that depends only on the underly- 
ing lattice. For parameters typical of polyacetylene this 
integral generates a factor of '^ 1.1, indicating a modest 
cooperative effect. However, the integral in Eq. [TT] is 
sensitive to the ratio A/2io and the cooperativity may 
be significant in other polymer systems. The transition 
from single site binding behavior to zipper-like binding 
occurs when the bound region grows comparable to the 
localization length ato/t' . 

Without the doping level to limit the length of the 
zipped region, the amount of chain in the tightly-bound 
and free states is determined by a competition between 
the decrease in the system's total electronic energy and 
the decrease in the chains' configurational entropy associ- 
ated with the binding regions. The statistical mechanics 
of the binding of two conjugated polymers is thus sim- 
ilar to that of DNA melting [26|. Here we adapt the 
Poland-Scheraga (PS) model designed for that problem 
to the current context. There are, however, significant 
differences between the two systems arising from the dif- 
ference in bonding mechanism and the lack of unique 
binding sites on the conjugated polymers, due to the ab- 
sence of the well-known DNA base pairing mechanism in 
the system of current interest. 

A pair of bound polymers consists of alternating re- 
gions of "zipped" (i.e. tightly bound) chains separated 
by loops of unbound polymer. The partition function for 
chains of length A^ in such a configuration may be written 
as 
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Here u(i}r , iJr ) and ^(jt) arc the Boltzmann weights of 
the tt"^ loop and the tightly bound domain respectively, 
enumerated from the left end of each chain. The number 
of sites in a tightly bound or "zipped" region, j^, must 
be an even number as each chain contributes j^r/2 sites. 



The size of the loop regions, i 



,•(1) , „-(2) 



are unconstrained 



because the translational symmetry of the chains allows 
the formation of asymmetric loops (i.e. the number of 

sites contributed by chain 1 , iJr , is not necessarily equal 
to the number of sites that chain 2 contributes to the 
loop, ijT )■ The sums over the number of "zipped" and 
unbound regions p and the length of each region are sub- 
ject to the constraint that the overall lengths of the chains 
are fixed 
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The statistical weight of a "zipped" region of n bound 
sites is given by v{n) = v"'a{n), where cr{n) accounts 
for the potentially length-dependent boundary energy be- 
tween looped and zipped regions, and 



is the Boltzmann weight associated with one bound pair 
of sites (one from each chain) in the "zipped" region. The 
energy that appears in the exponent has been broken into 
two parts. The first term E^ip is the per site inter-chain 
binding energy given by Eq. [Tljfor long tightly bound re- 
gions and — 2(4iQ -I- 1'^)^/^ 4- 4to for a single binding site. 
In practice we use the form given by Eq. [Tl]for "zipped" 
regions of length greater than ato/t' . Otherwise, we use 
the result for an isolated binding site. The second term 
represents a local per-site repulsion between the chains 
in close proximity. This energy arises from unfavorable 
steric interactions between the chains enforced by the lo- 
cal binding geometry and possibly electrostatic repulsion 
between the chains. The details of this energetic term 
will clearly vary with the chemical details of the spe- 
cific conjugated polymer system in question. Since Ees 
is of order ksT, we expect that the net binding energy 
may be smoothly varied from a minimum of Ezip to net 
repulsive values. The parallel configuration considered 
below may also contain a van der Waals component en- 
hancing the net attraction [13|. However, this interaction 
decays much more slowly with distance than interchain 
hopping, and therefore the details of aggregation on the 
monomer level will be dominated by the mechanism con- 
sidered here. We do not consider these details further 
here. 

The transition between a loop and a zipped re- 
gion of length n results in a boundary free energy 
— fc^rin (cr(n)). This boundary energy has contributions 
from the cooperativity of the intermolecular bonds, the 
configurational entropy cost associated with constraining 
the polymers' backbone with the formation of the sec- 
ond bond, and a bending energy associated with forming 
the "Y" -junction where the coUinear polymers split into 
an unbound loop. In the case of conjugated polymers, 
the first two contributions are negligible due to the weak 
cooperativity of the bonds and the limited fiexibility of 
the chains. This is in contrast to the situation in nucleic 
acids where the aromatic base stacking interactions leads 
to strong cooperativity in the inter-chain bonds while the 
flexible backbone results in a negligible bending energy 
and a non-negligible entropic contribution. 

For conjugated polymers the statistical weight for the 
bound regions takes the form 



v{n) =0 n = odd 
= V 11 = 2 
= cru" n > 4. 
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This form accounts for the fact that if the "zipped" region 
is only a single site long, i.e. a single binding site, the 
angle between the chains is only weakly constrained. For 
longer tightly bound regions, however, the chains must be 
parallel over its length then bend sharply at the start of 
the surrounding loops. The boundary parameter a may 
be estimated using the conjugation length in a melt. Un- 
der these conditions the polymer has conjugation breaks 
approximately every ten sites [27[. Since there must be 
two such breaks per tightly bound region and there are 



four ways of arranging the breaks, the overah Boltzmann 
weight should be ct ^ 0.04. 

The statistical weight for an unbound loop of the 
paired chains accounts for the configurational entropy of 
the self-avoiding loop when the loop's arc length is long 
compared to the thermal persistence length of the poly- 
mer. When the loop is short compared to this length, the 
bending energy of the chain makes the dominant contri- 
bution to the free energy of the loop. These two effects 
are captured using a weight of the form 
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with u = e ' *= . Here (.k is a length on the order of the 
persistence length ^^ ~ £p ^ a, and we have assumed 
that short loops trace out the arc of a circle [1^. The 
constant c accounts for the excluded volume of the chains 
and takes the value ^2.1 in three dimensions |28| . 

The sums in Eq. [T^] are difficult to evaluate because 
of the restriction imposed by Eq. [T31 They can be made 
more tractable by relaxing this constraint. In order to 
do so we study the function 
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which amounts to working in the grand canonical ensem- 
ble of the polymer length where we have introduced a 
fugacity x^^ . With the use of Eq. [12] this sum can be 
rewritten as 



r(x) = Lix)R{x)J2{Uix)V{x)Y 

p=0 
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where the functions U and V are defined by 
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The functions L{x) and R{x) are defined in analogy to 
Eqs. [T9]and account for the free energy of the tails at the 
ends of the chain. L{x) and R{x) will have a different 
functional form than the portions of chain on the inte- 
rior of the complex due to the fact that the chain may 
terminate in either a bound "zipper" or a pair of free 
tails that are not subject to a loop closure constraint. 
However, these details are unimportant as the free en- 
ergy of the system is dominated by the chain interior in 
the thermodynamic limit. 

In the limit of an infinitely long pair of chains i.e. 
where N — > cxd, the partition function Z{N) of the two- 
chain system must scale as x1^ , so that the free energy 
per site is now independent of chain length. The param- 
eter xi is as yet undetermined, but is related to the free 



energy per monomer / by / = —ksT \n{xi). From this 
we recognize the sum shown in Eq. [T7] will converge for 
all X > xi. We now determine a;i by examining conver- 
gence properties of this sum as given by Eq. 1181 Noting 
that the end effects L{x) and R{x) necessarily generate 
bounded corrections, we see that the divergence of the 
sum will occur at the roots of U{x)V{x) = 1; we will 
ignore all such end effects here and in the following. As 
we want the dominant term in the limit of long polymers, 
the free eiiergy is controlled by xi, the largest root of this 
equation [29|. From Eqs. 1151 and 116 1 we find that xi is 
then given by the solution of 
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This result allows us to determine the free energy and all 
thermodynamic properties of the paired chains. 

From this result we may compute any number of phys- 
ically measurable quantities involving the structure of 
paired conjugated polymers. We will consider three here. 
First, we compute the fraction of paired chains as a func- 
tion of polymer concentration in dilute solution, where we 
will show that, in the limit of sufficiently long polymers, 
there is a sharp cross-over between free chains and bound 
pairs as a function of concentration. One may imagine 
that in the strongly bound limit, more complex bundles 
of chains should form. Second, we examine the distribu- 
tion of free loops in the paired chains as a function of 
the strength of the binding interaction. This result may 
have measurable consequences for the diamagnetic sus- 
ceptibility of conjugated polymer solutions. Finally, we 
compute using a minimal set of assumptions the effec- 
tive persistence length of the paired chains as a function 
of the strength of the binding interaction. Such results 
have implications for the interpretation of scattering data 
from these polymers in solution. 

To calculate the fraction of bound chains we write the 
free energy for a dilute polymer solution 



VksT 



= ci (In (ci) — 1) — CiA^lnu 
+ C2(ln(c2) -l)-C2Nlnxi 



(21) 



in terms of the number densities of the unpaired and 
paired polymers, ci and C2 respectively. V is the total 
volume of the system. To find the fraction of bound 
chains as a funtion of the total polymer number density 
Co we use the conservation of the number of chains to 
write ci ~ cq ~ 2c2 and minimize Eq. [5T] with respect to 
C2. From this we find that 
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Because of the large polymerization index A^ ::^ 1, the 
transition from single chains to bound pairs occurs over 
a very narrow range of binding energies and becomes even 
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sharper as N increases [ij]. The binding curve is shown 
in Fig. ini 




FIG. 6: (color online) (top) The apparent persistence length 
as a function of interchain binding energy. For sufficiently 
small binding energy and for highly flexible (low persistence 
length) individual chains, the polymers are remain unbound 
and have an apparent persistence length £min equal to that 
of the individual chains. As the binding energy increases the 
polymers begin to form complexes, indicated here by the solid 
line which shows where exactly half of the chains are paired: 
2c2/co = 0.5. For strong enough binding energy, the average 
length of the zipped regions (L^ip) exceeds £max (dotted line), 
taken here to be 400a. Between the two lines the apparent 
persistence length is given by (Lzip) since the paired chains 
will be essentially straight in their zipped regions, but have 
free hinges associated with the looped ones, (bottom) Appar- 
ent persistence length plotted along the red line above (red). 
Also shown are curves for smaller and larger values of the 
nucleation parameter a. Although we show the persistence 
length saturating above Inu ~ 0.6, we expect that sufficiently 
strong binding will drive the formation of larger aggregates, 
further increasing the apparent polymer rigidity. 

The physical properties of bound chains will be de- 
termined in part by the statistics of the zipped and loop 
regions. For example, the magnetic susceptibility of poly- 
meric solution should depend in part on the number den- 
sity of semiconducting or conducting (for suitably doped 
chains) loops in solution, which itself will depend on the 
loop fraction of the individual bound chains. The frac- 
tion of zipped regions, which are significantly stiffer me- 



chanically, will determine the effective persistence length 
of the polymers. We first calculate the distribution of 
unbound loops. 

We consider a pair of bound polymers of length N in 
the limit TV ^' oo so that end effects may be ignored. 
The partition function for this system is 
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By translational symmetry all the loops are identical so 
we calculate the mth moment of the length of the first 
loop from 
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By the definition of xi, U{xi)V{xi) = 1, so we have 
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In Fig [7] we plot the average loop size along with its 
variance. The loops are typically just a few persistence 
lengths as this length is required to avoid a prohibitive 
bending energy penalty. The narrow distribution is in- 
dicates that longer length loops also incur a prohibitive 
binding energy penalty. 




FIG. 7: The average loop length {Lioop} (solid line) is plot- 
ted as a function of the per site binding strength of the 
chains using a bare persistence length of the noninteract- 
ing chains of £p = 10a. The dotted lines provide a mea- 
sure of the width of the distribution of loop sizes by plotting 



(Lloop) ± \J {Lfoop) ^ {Lloop)'^- 



We now turn to calculation of the effective persistence 
length of the paired chains. The local bending modulus of 
the polymer changes markedly when to chains are tightly 
bound. One should thus be able to distinguish free loops 
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from "zipped" regions of paired chains purely from their 
local conformational statistics in equilibrium. We first 
demonstrate this difference in persistence length by esti- 
mating its value for free chains and tightly bound ones. 
In their free state the persistence length of the polymers 
may be inferred from the conjugation length observed in 
polymer melts. Given a finite conjugation length of this 
form, one may simply describe the conformational statis- 
tics of the free polymer as a freely rotating chain having 
a segment length equal to its conjugation length. This 
length imin is on the order of ten monomers at for the 



Eq. 



We obtain 



case of PPV [27[ , but this value depends on the specific 
chemical system. 

When the chains are tightly bound together in what 
we call the "zipper" state, the intermolecular bonds will 
prevent rotations about the polymeric backbone. We 
imagine the principal remaining source of backbone flex- 
ibility on the microscale to come from the bending and 
stretching of the backbone chemical bonds. The en- 
ergy scales associated with bond deformation are sig- 
nificantly higher than for the bond rotation of the free 
chains. The bond stretching modes of the backbone are, 
in fact, too stiff to play a significant role. Using the es- 
timation of the effective bond-stretching modulus of 21 
eV/A^ we estimate that the effective paired-chain mod- 
ulus for bending perpendicular to the conjugation plane 
to be ~ 120 eV/radian [2|. This is to be compared with 
the 5.5 eV/radian modulus for bending the bonds in the 
conjugation plane [301 • We estimate that effective per- 
sistence length of the paired chain by determining length 
over which a one radian uniform bend generates one ksT 
of clastic energy and find that the paired-chain persis- 
tence length is at least one order of magnitude bigger 
than that of the free chain, i.e. £max — lO^min- 

A representative equilibrium configuration of the 
paired chains will, in general, consist of a series of alter- 
nating "zipped" regions and loops in which the two chains 
are locally unbound. Since the "zipped" regions are sig- 
nificantly stiffer than the floppy loops, we may approxi- 
mate the resulting configuration of the paired chains as 
a random walk consisting of essentially straight "zipped" 
regions connected by flexible loop regions at which the 
polymer bends freely. To determine the effective persis- 
tence length of the paired chains then we must determine 
the mean length of the stiffer "zipped" regions in ther- 
mal equilibrium. This length is given by the ratio of frac- 
tion of bound sites 9 to the fraction of sites that make 
up the boundary of a "zipped" region: Oend- The latter 
value is calculated by appending a boundary term to the 
loop regions such that U{x) —^ auU{x) that counts such 
boundaries. The fraction of sites that incur this bound- 
ary weight is then given by 
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which may be evaluated by implicitly differentiating 
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The parameter 9, determining the number of bound sites 
is given by the derivative 
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Alternatively, 9 could be calculated from an expression 
analogous to Eq. [5S1 The average length of the "zipped" 
regions is 
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The various regimes of the effective persistence length 
can be seen in Fig. |6] as a fraction of the persistence 
lengths irainAmax and strength of the binding interac- 
tion as parameterized by Inti - see Eq. [TD Based on 
our calculations, we expect this last parameter to be of 
order one. In this figure we see that the effective per- 
sistence length smoothly crosses over from Irnin to ^max 
as the binding interaction increases for sufficiently long 
polymers, i.e. for those polymers whose contour length 
L is greater than (.max- 



VI. DISCUSSION 

We have further explored the effect of inter-chain elec- 
tronic tunneling on polymer pairing (i.e. forming two 
chain bound states) and on the structure of those pairs in 
thermal equilibrium. At the level of the formation of an 
isolated binding site, the binding energy is highly insen- 
sitive to the density of states at the Fermi energy of the 
respective chains. In particular, doped metallic chains, 
explored earlier [13|, and semiconducting ones form in- 
dividual binding sites that decrease the total electronic 
energy of the system by O(t'^/to) ~ ksT. We also inves- 
tigated the electronic interaction between tunneling sites 
that are separated by a short arc length along the chain. 
Here we see that there is a weak attractive interaction 
between binding sites due to the interaction of the lo- 
calized states at these tunneling sites. The strength of 
this attraction is, for at least one physical set of model 
parameters |12l |. rather weak, ^ O.lksT, but this value 
may vary widely between conjugated polymer systems. 
The attractive interaction may in some cases lead to the 
clustering of bound regions. 

We support the single tunneling site results of our 
simple tight-binding model using Hartree-Fock numerical 
calculations; the numerics also show that the inter-chain 
binding is somewhat insensitive to the precise orienta- 
tion of the chains at the tunneling site, although the an- 
gle between the chains at the crossing point does lead to 
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significant changes in the binding energy. There we find 
that the perpendicularly crossed-state which simultane- 
ously allows for the a small separation at the tunneling 
site and minimizes inter-chain steric repulsion is actu- 
ally not the lowest energy state. This configuration is 
a local minimum, while the global minimum occurs for 
more nearly parallel chains in an orientation consistent 
with the alternating binding and non-binding sites corre- 
sponding to the ground state of the tight-binding model. 
In the context of that model, we understand the forma- 
tion of the "every-other" structured ground state in terms 
of enhancing the Peierls splitting of the conduction and 
valence bands. It appears from the Hartree-Fock calcula- 
tions that the local structure of a single binding site also 
favors the formation of such a "every-other" bound state. 
It is interesting to note that the crystal structure of poly- 
acetylene shows evidence of this "every-other" bonding 
pattern |3ll. |32|. In this case, however, each polymer in 
the crystal makes contact with several others, but the 
local equivalence of the pattern of contacts suggests that 
the electronic mechanisms we consider contribute to the 
stabilization of this crystal structure. 

The binding between two chains is enhanced by the 
co-localization of solitons with the tunneling sites. The 
mixing of the mid-gap states associated with the solitons 
by the tunneling matrix element leads to this binding 
enhancement for uncharged solitons, but not for charged 
ones. This effect is weak when there is only one soliton 
at a tunneling site. Thus, solitonic binding enhancement 
should scale as the square of the equilibrium (uncharged) 
soliton density. Recalling that their density in thermal 
equilibrium is low, we do not expect this to play a large 
role. It is important to note, however, that a binding 
event between a charged soliton and an uncharged soli- 
ton will result in one of the gap states being half filled 
and therefore lead to half of the total possible binding en- 
hancement. The density of charged solitons is easily con- 
trolled by doping. Using this, one could, in the presence 
of a sufficient density of charged solitons, produce a sig- 
nificant solitonic enhancement (on the order of Ifc^T per 
site) of the chain binding energy that scales linearly with 
the density of uncharged solitons. We expect the same 
behavior to be observable in polaron/tunneling site inter- 
actions in conjugated polymers more chemically complex 
than polyacetylene. 

The statistical mechanics of two paired conjugated 
polymers is rather subtle since the binding sites generate 
free energy changes of only ~ Ifc^T per site and allow 
for significant local conformational freedom at a binding 
site. The structure of the paired chains results from a 
nearly equal competition of chain configuration entropy 
and binding energy. Of course, many of the same sub- 
tleties have been addressed in the problem of DNA melt- 
ing. The one main difference between DNA melting and 
the binding of conjugated polymers, is that in the lat- 
ter there is no equivalent of the base-pairing mechanism 
that promotes the binding in registry of complementary 
strands. The lack of pairing registry in the current prob- 



lem enhances the entropic gain of loop-formation by al- 
lowing loops of a fixed length L to be created from various 
lengths li and I2 (such that h + I2 = L) from the two 
chains. Using a modified Poland-Scheraga model which 
accounts for this difference, we found that the effective 
persistence length varies over at least an order of magni- 
tude as a result of inter-chain binding. The bound regions 
are significantly stiffer than the unpaired or free loops of 
the chains allowing us to calculate the effective persis- 
tence length of the paired chain in terms of the mean 
length of the stiff, bound or "zipped" regions. 

A number of open questions remain. The spectroscopic 
signature of the pairing mechanism on two chains is still 
unresolved. Moreover, the formation of polymeric ag- 
gregates of more than two chains remains an open ques- 
tion, particularly with regard to the formation of crystals 
of conjugated polymers. The fundamental tight-binding 
approach to the electronic degrees of freedom coupled to 
the statistical mechanics of the chains, further developed 
in this article, should provide the basis for these future 
investigations. 
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Appendix A: Transfer Matrix Calculations 

The one dimensional character of the electronic hop- 
ping mechanism lends itself naturally to a transfer ma- 
trix approach. In this appendix we review the formal- 
ism as it is used to solve for the electronic eigenstates 
of interacting polymers. In the unperturbed state the 
eigenstates are normal Bloch waves. The introduction of 
the inter-chain interaction will scatter these waves and 
lead to localized states. We begin from the electronic 
part of the interacting Hamiltonian, Eqs. [1] and [3l and 
work in the basis of chain symmetrized states as shown 
in Eq. 2] After this transformation, the symmetric and 
anti-symmetric parts of the Hamiltonian may be written 

Hs/A = -^i£/+i(K+l,a)(^,aKK,a)(^+l,a|) 

+A|0,a)(0,a|, (Al) 

where A = — t'(+t') in the symmetric (anti-symmetric) 
sub-space. 

We solve for the energy eigenvalues E and the eigen- 
vectors \E) = ^^ ci\t) that satisfy the Schrodingcr equa- 
tion Hs/a\E) = E\E). By operating on the Schrodingcr 
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equation with the vector (€| we find that the amphtudcs 
ci obey the set of simuhaneous equations 



— ti-ijC£-i + £lCl — t£j+iC£+i = EC£, 



(A2) 



which, along with a trivial identity, may be written in 
matrix form. 
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From Eq. lAll the on-site energies ei vanish everywhere 
except at the binding site, e^ = SqiX. 

To bring the Hamiltonian into Bloch form, we con- 
struct the transfer matrix for a unit cell 
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If we number the carbons according to 1=2-3=4-5. . . , 
where the single and double lines represent long and 
short bonds respectively. Eg. I A4I propagates of the wave- 
function between consecutive odd sites. As required by 
the discrete translational invariancc of the Hamiltonian, 
Eq. IA4I takes a diagonal form in the momentum basis 
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where 2a is the unit cell dimension of the dimerized chain, 
and 
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The first of these equations implicitly redetermines the 
dispersion relation of the Bloch waves making up the en- 
ergy eigenstates of the system. The solution E{k) is given 
below in Eq. IA91 For completeness we also report the 
change of basis matrix defined by Eg. lASl in terms of the 
dimensionless parameter q given by 



o'9 



-E 



l^g-ilka j^l^ 



(A8) 



The complex number is the phase change of the wavefunc- 
tion associated with the quantum number k (wavenum- 
ber) across the long bonds separating one unit cell from 
the next. In other words: q = Q+ie"', as may be checked 
using Eas. lA31IA41 and lASI The phase shift within a unit 
cell is 2ka — g so that the total phase shift from one unit 
cell to the next is 2fca, as required by Bloch's theorem. 

Of course, the condition for bounded eigenfunctions is 
that k must be real, leading to the dispersion relation 



E{k) = ± 



tl+tl 



2tit2 cos{2ka) 



(A9) 



which has a total bandwidth of 4to and a bandgap cen- 
tered around E = oi width 8uoa = 2A. 

To study the effect of the binding site on the electronic 
spectrum we repeat the transformation Eq. IA5l with the 
dimer containing the binding site at £ = 0. The result is 



Mx = UMiMqU' 



i2ka 



e'-^'il+iS) 



A2ka 



— iSi 



, — i2ka 



*2fe'^(l-i(5) 

(AlO) 
where S = XE /2tit2 sin{2ka). In addition to the states 
given by Eq. [5] for — tt < 2ka < tt, the off-diagonal ele- 
ments in Eq. lAlOl allow normalizable states with complex 
wavenumber through the scattering of a growing wave 
into a decaying one. This is accomplished when the up- 
per left matrix element vanishes, therefore 



1 



-i5 



-iXEh 



2tit2 sin(2fca) 
XEt 



ViE't-tf-ti) 



W 



~Wi' 



(All) 
(A12) 

(A13) 



With a little more algebra wc find that the four bound 
states shown in Fig. [1]: are given by 



^'±A/4f2t2+A2(t2 + f2)+^' 



Eb = ±sm^{X)\l tl+tl , 2 -y-i'-2 ' - vn ' -2/ ' 4' 

(A14) 
where the upper sign applies to the "ultraband" states 
(Eq. [S]) and the lower sign applies to the gap states 

(Eq.Ei im. 

We now compute the analogous bound state as- 
sociated with a tunneling site localized at the cen- 
ter of an idealized soliton. We represent the soliton 
as a single-site disruption of the dimerization pattern 
as discussed in Section IIVI It may be written as: 
o = o — o = o — • — o = o — o = o and has undis- 
torted unit cells on either side of the soliton's center (•) 
connected to that special site by the long bonds (— ) that 
have tunneling matrix element ^2- The tunnel matrix 
element also occurs at this special and central site. 

The transfer matrix at the central site is 



M, =\ X-E 

t. 



Upon changing basis as shown in Eq. IA5I wc find 
1 



UALU-^ 



2i sin(fc) 

2 ^~iQ \-E I ^ p-2iq 



(A15) 



(A16) 



'JQ \-E 



iq \ — E o I iq X—E 

t2 t2 



Once again, the bound states exist for values of E where 
the upper left element vanishes, so we have 



1 



-iq 



X-Eb 

2t2 



(A17) 
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Using Eqs. |Ml ESl we find 

E^ — XEb 



1 



(A18) 



E^-q + t^, ^{Ei-ii-iiy-Aiiii 

The solutions are given by the roots of the polynomial 

i) = Et^ El{Ul + A2 + a2) + 4£;bAtoA. (A19) 

It is easily verified that the solution _E(, — does not 
satisfy Eq. IA17[ and thus we arrive at Eq. [H 

Appendix B: Hartree-Fock Calculations 

To better understand the limitations of the tight- 
binding model and explore the role of steric interactions 
near the tunneling site for the specific case of polyacety- 
lene we numerically computed the binding energies of two 
short (20 monomers each) polyacetylene oligomers with 
quantum chemical calculations using the package Gaus- 
sian [33|- After relaxing a single chain of fixed length, 
we introduced an identical chain in various geometries. 
In each case the interaction energy was calculated as the 
energy of the two chain system minus twice the energy 
of the single chain. Energies were evaluated using the 
non-counterpoised Hartree-Fock technique with the 6- 
31G* basis set. Additional calculations using the B3LYP 
method found results qualitatively similar to those from 
Hartree-Fock (data not shown). 

Using this procedure we tested the distance and orien- 
tation dependence of the binding energy due to a single 
tunneling site. We parameterized the geometry of the 
near miss of two chains by two angles and one distance. 
As shown in Fig. [51 the perpendicular distance between 
the two chains at the binding site is h. The angle (j) mea- 
sures the rotation of the plane containing the intra-chain 
carbon bonds of one polymer about an axes parallel to 
that chain. Finally, the angle between the long axes of 
the two chains (looking down along the direction parallel 
to the vertical displacement h) is given by 6. In Fig. [9] 
we show the binding energy obtained numerically as a 
function of h and 6 for a fixed value oi (f) — Q (i.e. with 
the conjugation planes parallel). As expected for the 
tunneling mechanism, the attractive interaction decays 
rapidly (exponentially) with distance reaching minimum 
at h = 4.0 A for all values of 6. The angular dependence 
at fixed h is more complex. For h = 4.0A correspond- 
ing to maximum binding, the attractive interaction has 
a minimum at the perpendicular crossing of the chains 
{6 = 7r/2) where the net binding energy is —0.0137 eV. 
This local minimum is globally unstable; there are deeper 
minima near 9 — 7r/6, 57r/6 corresponding to more nearly 
parallel alignments of the chains. The 9 ~ 5tt/6 min- 
imum is deepest and in closer agreement to the "every 
other" site ground state predicted by the tight binding 
model. For precisely parallel alignments {9 = 0) in which 
the chains are stacked on top of each other as shown in 





FIG. 8: (top) Two polyacetylene chains from the view (top) 
showing the vertical separation between the chains h and axial 
rotation angle </>. The tunneling site is shown as the (green) 
pillar, (bottom) The same chains from the top view, looking 
down along the direction of h and demonstrating the crossing 
angle d. 



Fig. El they do not attract, as predicted by our tight bind- 
ing model. In fact, interactions not included in our model 
give a repulsive interaction of 0.154 eV. Our results are 
consistent with previous work that has shown the maxi- 
mally overlapped configuration to be a high energy state 
with more favorable interactions occurring when one of 
the molecules is rotated or translated such that the over- 
lap is removed from some of the carbons [33, ISa] . 

The binding energy shown in Fig. [5] is also insensitive 
to the precise lateral alignment of the two chains at the 
binding site. We shifted the top chain by 0.2 A in the 
positive and negative X and Y directions, where h lies 
along the Z direction. These shifts altered the binding 
energy by less than 6%. These data are not shown. 

Finally, we rotated the upper chain by the angle cf) as 
defined in Fig. [8l This rotation preserved the 4 A sepa- 
ration between the two carbons at the binding site. For 
negative rotations the interaction becomes less favorable 
due to the steric repulsion of the nearby carbons of the 
two chains, but for positive rotations the binding energy 
increases to 0.0144 eV at an angle of = 0.3, as shown 
in Fig. Uni Note that these energy shifts are significantly 
less than fc^T suggesting that the rotational motion of 
the polymers is unconstrained by the tunneling sites in 
thermal equilibrium. 

We find that these results arc only weakly dependent 



15 




150 



FIG. 9: Binding energy between two 20 carbon polyacetylene 
chains as a function of their separation h and contact angle 
6. The optimal separation is h — 4.0A at all angles. There 
are pronounced local minima at 6 — 30° and 6 = 90° and 
a global minimum at 6 = 150°. The flat surface represents 
regions where the net binding energy is repulsive within the 
HF approximation. 
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FIG. 10: The interaction energy as a function of </> for h = 4.0 
A and 61 = 90°. 



upon chain length; for example, as the chain lengths were 
increased from ten to forty carbons, the depth of the 
well increased by about five percent. In order for the 
numerically obtained binding energies to match our tight 
binding results, we need to set to ~ 2.5 eV and t' ~ 0.26 
eV. This result is consistent with the overlap measured 
in Ref. [12| for configurations where the contact extends 
over the entire molecule. 

The Hartree-Fock calculations presented here support 
the notion of intermolecular bonds forming at the close 
approach of two polyacetylene chains. The strength of 
these interactions is consistent with previous studies J12l | . 
The range of the attraction is as expected for the pro- 



posed tunneling mechanism and its angular (9) depen- 
dence suggests both perpendicular and nearly parallel 
arrangements are favored for a single binding site. The 
binding interaction, however, is sufficiently insensitive to 
the precise orientation of the polymers as to allow the 
polymers to form a variety of aggregates. 

We conclude by noting that the ab initio calculation 
of the interactions between 7r-conjugated molecules is an 
area of ongoing research [3J, [Sa] . While significant ad- 
vances have been made for small aromatic molecules, 
these molecules lack the extended orbitals essential to 
our model. These methods are not feasible for the larger 
systems we require, and therefore, we have employed the 
significantly less accurate HF method. While we are en- 
couraged by the qualitative agreement with our tight- 
binding model, we emphasize that our HF results should 
not be viewed as quantitative. 



Appendix C: Binding Site Interactions 

To examine the interaction energy between two tun- 
neling sites we consider a chain of A^ tight binding sites. 
We place the two tunneling sites, which act like impurity 
potentials, at sites M and M + d such that the the tun- 
neling sites are symmetrically placed about the center of 
the chain: N = 2M + d — I. This provides a computa- 
tional convenience without greatly reducing the validity 
of the results so obtained. The position of the localized 
states along the chain should be irrelevant on long poly- 
mers A^ ^ 1 provided neither tunneling site is within its 
localization length (a few lattice constants) of the chain 
ends. 

Since the even-odd effect observed in Fig. [2] is un- 
changed by the chain dimerization, we consider the case 
of a uniform chain with i^^^+i = io and dispersion 
E{k) — —2tQC0s{k). The eigenstates have the form 



A/-1 



|fc) 



y^ sin(H)|£) 

£ 1 




M+d 




+A J2 sm{ke + (i))\e) 




i=L 




N 




± Y, sin(A:(iV + 1 - 


-eW) 


e=M+d+i 





(Cl) 



where A is an amplitude to be determined, and the sign 
of the final term is determined by parity. The boundary 
conditions require that the wavefunction be an eigenstate 
of the parity operator and and remain continuous at the 
impurities so that 



kd 
2 



2 



kM 

sin(fcA/) = ^sin(fcM 4- 0). 



(C2) 
(C3) 
(C4) 
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We derive a final condition by requiring that the 
Schrodinger equation is satisfied at the impurity sites 

-2toCOs(fc)sin(fcM) = Asin(fcA/) 

-iosin[A:(M - 1)] - taAsm[k{M + !) + (/)]. (C5) 

Using Eqs. IC2| IC3[ and lCSI wc find the conditions on k 
for states of even and odd parity are 

A kd 

--— = tan( — ) — cot(A:M) even states(C6) 

to sin(fc) 2 

-— = - cot( — ) - cot(kM) odd states. (C7) 

iosin(A:) ^ 2 ' ^ ' ^ ' 

We would like to examine the change in energy of the 
states as the distance between the impurity potential is 
changed but the length of the chain is constant. To do 
this we make the substitution M = {N — d + l)/2 and 
find that the quantization conditions become 



cot k 



tan k 



iV+l\ 


^[l + cos{kd)] 


2 ) 

N+l\ 


sin(fc) — 27" sm(kd) 
^[l-cosikd)] 


2 ) 


sin(fc) + 27;^ sm{kd) 



even(C8) 



odd.(C9) 



To see the energy shift of the individual states it is help- 
ful to graphically solve these equations. In Fig. [Tl] we 
plot both sides of Eq. IC8I for the symmetric and anti- 
symmetric states (A/2to = ±0.5) and d = 8. The left 
side of the equation, drawn in dashed lines, is a series of 
tangent curves with a spacing set by the length of the 
chain. For very long chains these tangent curves will be 
spaced infinitesimally close together. The allowed k val- 
ues are found at the intersections of the solid curves and 
the dashed curves. In the limit A — !■ 0, the solid curves 
lie along the k axis and we find the unperturbed states 
satisfy k = (2n -I- l)7r/2(7V -I- 1) where n is an integer. 
For small positive values of A the right side of the equa- 
tion takes positive values and the eigenstate is shifted to 
larger k values. Similarly, negative values of A shift the 



states to longer wavelengths. Therefore, the amplitude 
of the solid curves in Fig. [11] grows with the magnitude 
of the wavenumber shift of the eigenstates. 

We see in Fig.[TT]that, although the curves for positive 
and negative values of A oscillate with the same period, 
the curves are not symmetric about the k axis. Starting 
at small fc, there are alternating regions where the sym- 
metric states are shifted more than the anti-symmetric 
states and vice versa. The origin of this asymmetry be- 
tween the sub-spaces is the phase shift induced by the 
impurity potentials discussed in section fllBl The contri- 
bution of these states to the total energy of the system 
will depend on the location of the Fermi level. If the sep- 
aration between the binding sites, d/a, is even the Fermi 
level, Ef = 0, lies just above a region of large shift for the 
symmetric states and below the corresponding region of 
large shifts for the anti-symmetric states. This results in 




FIG. 11: Graphical solution of Eq. ICSK see text). 



a net negative contribution to the binding energy. How- 
ever, if the separation is odd the Fermi level lies immedi- 
ately above a region of large shifts for the anti-symmetric 
states and the contribution is positive. 
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